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According to Shoemaker, the "impact of solid bodies is the most fundamental process that has 
taken place on the terrestrial planets" [l|, as they shape the surfaces of all solar system bodies. A 
lot of information on this process has been extracted from remote observations of impact craters 
on planetary surfaces. However, the nature of the geophysical impact events is that they are non- 
reproducible. Moreover, their scale is enormous and direct observations are not possible. Therefore, 
we choose an alternate and of course downscaled experimental approach in order to guarantee 
reproducible results: We prepare very fine sand in a well defined and fully decompactified state by 
letting gas bubble through it. After turning off the gas stream, we let a steel ball fall on the sand. 
The series of events in the experiments and corresponding discrete particle simulations is as follows: 
On impact of the ball, sand is blown away in all directions ("splash") and an impact crater forms. 
When this cavity collapses, a granular jet |3,l3l emerges and is driven straight into the air. A second 
jet goes downwards into the air bubble entrained during the process, thus pushing surface material 
deep into the ground. The air bubble rises slowly towards the surface, causing a granular eruption. 
In addition to the experiments and the discrete particle simulations we present a simple continuum 
theory to account for the void collapse leading to the formation of the upward and downward jets. 
We show that the phenomenon is robust and even works for oblique impacts: the upward jet is then 
shooting backwards, in the direction where the projectile came from. 



It has long been known that jets can be created when 
a ball or a fluid droplet impacts on a fluid surface 
d HE nil- S. Thoroddsen and A. Shen found similar 
jets on impact of lead spheres on monodisperse spheri- 
cal glass beads [3- We did similar experiments on fine 
sand, but found it hard to achieve quantitatively re- 
producible results, presumably due to the random na- 
ture of the force-chain-networks in the granular material 
13, [l3, HB US' Therefore, in order to prepare a well- 
defined initial state, we decompactify and homogenize 
extremely fine sand (average grain-size of about 40/im; 
grains are non-spherical) by blowing air through it via 
a perforated bottom plate. The height of the sand bed 
above the bottom plate is typically 25-40cm. The air is 
slowly turned off before the experiments and the grains 
are left to settle in an extremely loose packing with the 
force-chains either broken or substantially weakened. We 
call this a "fluid-like" state. Impact events on this well- 
prepared fine sand will be gravity-dominated. We let a 
steel ball (radius i?o = 1.25cm) fall from various heights 
(up to 1.5m) onto the sand and observe the dynamics of 
the sand with a digital high-speed camera (up to 2000 
frames per second). 

The series of visible events is as follows (see figure^: 
First, the ball vanishes in the sand and a crown-like splash 
is created. Inhomogeneities develop in the crown, due to 
the inelastic particle-particle interaction (figure 1, frames 
3-5). Then, after a while, a jet shoots out of the sand 
at the position of impact. In all our experiments the 
jet height exceeds the release height of the ball. (The 
jets of ref. Q never reach the release height, because the 
sand is less fine and much less decompactified.) While 
the upper part of the jet is still going upwards, in the 



lower parts the inelastic particle-particle collisions lead 
to density inhomogeneities in the jet (figure 1, frames 7- 
8) . These inhomogeneities resemble those of the surface 
tension driven Rayleigh-instability of a water jet, even 
though there is no surface tension in granular matter. 
Finally, after about half a second, a granular eruption 
is seen at the position of impact, resembling a volcano 
(figure 1, frames 8-9). The collapsing jet first leaves a 
central peak in the craterjS^, but the granular eruption 
violently erases this peak. 

How does the jet form? To find out what is going 
on below the surface of the sand, we (i) performed di- 
rect numerical simulations, (ii) redid the experiments in 
two dimensions, meaning that we replaced the ball by a 
cylinder (with axis parallel to the surface and orthogonal 
to the side plates) which we let fall into a bed of sand 
between two transparent plates, and (iii) employed the 
analogy to jet formation in fiuids [1, IE H, 0, IE Il3 • 

(i) In the discrete particle simulations, the sand parti- 
cles are modeled as spheres which interact via inelastic 
"soft-sphere" collision rules. The interaction of the par- 
ticles with the surrounding air is included via empirical 
drag force relations |lq|. Since the maximum number of 
particles that can be simulated is presently of the order of 
one million, we can perform only quasi-two-dimensional 
simulations, where the thickness of the sand bed between 
the parallel plates is eight grains i33]. Altogether, the cal- 
culation includes N= 1.3 million homogeneous beads of 
density 1000kg /m'^ and diameter 500^m (i.e., approxi- 
mately a factor 10 larger than in experiment) in a con- 
tainer of 24cm X 0.4cm ground area and a sand bed 
height of about 17cm. The beads are pre-fluidized with 
air, just as in the experiments, and then a 1.5cm diameter 
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FIG. 1: Jet formation after the impact (vq = 2.43m/s) of a 
steel ball of _Ro = 1.25cm on loose very fine sand. The jet in 
this experiment exceeds the release height of the ball. Frames 
2-4: splash; frames 5-6: a jet emerges; frame 7: clustering 
within the jet; frames 8-9: granular eruption at the surface. 
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ball of density 3500kg /nt^ is dropped onto the beads with 
an impact velocity of 2 ni/s. The series of events can be 
seen in figure |21 revealing the jet formation process invis- 
ible in figure ^ The impacting ball creates a void which 
is then pressed together through the "hydrostatic" pres- 
sure from the side. [2J| At small depth the ball passes 
early, meaning an early start of the void collapse, which 
however is weak due to the small "hydrostatic" pressure. 
Conversely, at larger depth the collapse of the void be- 
gins later, but is stronger due to the larger "hydrostatic" 
pressure. Somewhere in the middle the collapse is fin- 
ished first, and the void walls hit each other. It is this 
singularity which leads to the formation of two jets: One 
upwards and one downwards into an air bubble which was 
entrained in the sand by the void collapse. The falling jet 
often leaves a central peak in the crater (which in our 3D 
experiments with the fine, decompactified sand is sub- 
sequently erased again by the granular eruption). Note 
that the jet in the discrete particle simulations is much 
less pronounced than in experiment. First, because the 
beads in the simulations are much larger than the sand 
grains in the experiment, i.e. the sand bed is less fluid- 
like and allows for less fine structure. Second, the sin- 



FIG. 2: Cut through the quasi-2D discrete particle simulation. 
Frames 1-3: the impact of the disk on the particles; Frames 4- 
6: the collapse of the void; Frames 7-8: the upward jet (which 
is less pronounced than in the 3D experiments). 



gularity due to the focussing along the axis of symmetry 
is weaker in 2D and quasi-2D experiments or simulations 
than in 3D, and the jet takes the form of a sheet. 

(ii) We performed such 2D jet formation experiments, 
by letting a cylinder fall into decompactified sand be- 
tween two transparent plates, and observing the jet for- 
mation process from the side (see supplementary mate- 
rial and ref. 3]). These experiments confirm the above 
sketched series of events. Again, the jet is less pro- 
nounced than in the 3D experiments. The entrained air 
bubble slowly rises in these experiments, finally leading 
to a granular eruption at the surface, just as observed in 
3D. 

(iii) The same series of events is also found after an 
analogous impact of a steel ball or a falling disk on water 
milallllllllll- We will employ this analogy below 
in order to set up a theoretical model. 

Before we do so, we discuss the role of the ambient 
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FIG. 4: Sketch of the void coUapse. When the accelerated 
sand grains from the sidewaUs of the cyhndrical cavity colhde 
on the axis of the cavity, two jets are formed: One downward 
into the entrained air bubble formed above the sphere, and 
one upward straight into the air. 
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FIG. 3: Cross-section of the 3D- void collapse following from 
our Rayleigh-type model, for the same impact velocity and 
ball radius as in figure 1. The void is pressed together by the 
"hydrostatic" pressure from the side, leading to a singularity 
and an upward and downward jet. 



air for the jet formation. We redid the discrete particle 
simulations with an air pressure reduced to nearly zero 
(vacuum), giving nearly indistinguishable results for the 
jet formation. The ambient air, however, can play a role 
during the evolution of the jet, provided that the impact 
velocity vq is very high. For (3D) experiments with very 
high impact velocities we observed that after the splash 
the crown goes inwards rather than outwards, due to the 
pressure reduction behind the fast projectile (Bernoulli's 
law). The crown in fact can fully close and the jet then 
hits the closed crown, leading to an explosion-like colli- 
sion (see the supplementary material) which spreads ma- 
terial all over the place. 

To work out the essentials of the void collapse, we 
now construct a "minimal" continuum mechanical model. 
First, the delay curve z(t) of the ball in the sand can 
be obtained from a simple force balance model involving 
drag, gravity, and added mass. It describes the experi- 
mental results obtained for a falling ball equipped with a 
thin tail rod, which allows for easy depth measurements 
[l9.] . The delay curve z{t) of the ball is inverted to obtain 



tpass (z), the time when the ball passes the layer of sand at 
depth z. This sets the initial conditions for the collapse 
of the two-dimensional void, namely R{z, tpass) = Ro 
and R{z, tpass) — 0. Here, i?(z, t) is the time and depth 
dependent radius of the void, see figure 21 

Next, the collapse of the void formed by the ball has 
to be described. It is driven by the ("hydrostatic") sand 
pressure p(z) at depth z. For small z the pressure simply 
is p{z) = Psgz, for larger z it saturates [13 . Here, ps is 
the sand density, assumed to be constant. If we neglect 
the dissipative processes both between the different layers 
of sand and between the sand grains in one layer, the 
dynamics for fixed depth z is determined by the Euler 
equation, 

Ps{dtv{r,t) + v{r,t)drv{r,t)) = -drp{r,t). (1) 

Here, v{r, t) is the velocity field in the sand. With conti- 
nuity dr{rv{r,t)) = 0, and with the boundary conditions 
v{R{t),t) = R{t) at the void's wall and v{Roo,t) = far 
away from the void, one obtains a Rayleigh-type [23, UM 
ordinary differential equation for each R{z, t), namely 



{RR + R^)log^ + l-R^^-p{z) 

-Koo ^ Ps 



gz. 



(2) 



The radius Roo is of the order of the system size, but the 
results only weakly (logarithmically) depend on this pa- 
rameter. The dynamics following from this Rayleigh-type 
model is shown in figure 13 resembling the void collapse 




one obtains that the time of the collapse t, 
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FIG. 5; Impact of the steel ball on soft, loose sand under 
an angle of approximately 45 degrees. Frames 2-4: forwardly 
directed splash; frames 4-5: a backward jet emerges; frame 6: 
clustering within the jet. 



in the discrete particle simulations figure [3 in the 2D 
experiments (see the supplementary material), in exper- 
imental work on the void collapse in transparent fluids 
13, H, H, Q, H, U^ , in boundary integral simulations of the 
complete hydrodynamical equations il3|i ^nd therefore 
presumably also in the 3D experiments in sand shown in 
figure n Just before and at the singularity {R{t) = 
and diverging velocity), the dynamics is determined by 
RR + R^ = 0, which has the solution R{t) ~ [ts - tf^, 
where t^ is the time of the singularity. The velocity there- 
fore has a square-root divergence R{t) ~ (ig — t)~^/'^. 

Having shown that the void collapse is driven by "hy- 
drostatic" pressure, we now can deduce scaling argu- 
ments [22 , for the limiting case of large impact velocity 
fo, which is the relevant one in the geophysical context. 
The time up to void collapse in depth z is the sum of 
the time z/vq it takes the ball to get there and the col- 
lapse time itself, which scales as ~ RqJ ^fgz. The depth 
Zc where the walls of the void first touch (i.e., the posi- 
tion of the singularity) can be obtained from minimizing 
this sum with respect to 2;, resulting in Zc/Ro ^ Fr^/^, 
where Fr = VQ/[gRf)) is the Froude number. From this 



tc - {Ro/vo)Fr'l^ ^ y^R^Fr-^/^ [2 
these scaling laws are consistent both with our contin- 
uum model and with our discrete particle simulations. 

We now come to the question how things change under 
an oblique impact J35j? We performed experiments under 
an angle of 45°. All other experimental conditions are as 
in figure ^ The series of events can be seen in figure 
Remarkably - but consistent with our theoretical model 
- the jet now points backwards, along the void created by 
the impacting ball. The backwards jet is also observed 
in the discrete particle simulations (see the supplemen- 
tary material). Note that this is different for an oblique 
impact on water, where the jet still goes upwards. 

We conclude the paper with speculations on possible 
implications of our findings on the im pac t mechanism 
within the geophysical context [11 IH lH |25;, 26, 2';^. 
However, we would like to caution the reader because - as 
pointed out above - lack of reproducibility of the details 
typify geophysical events. Moreover, though our decom- 
pactifying procedure minimized the relative energy stored 
in the ground as compared to the energy of the projec- 
tile, the absolute energy scales in our experiments are of 
course very different as compared to geophysical events. 
Nonetheless, we believe that the following speculations 
may stimulate discussions in a geophysical context: (i) 
After the impact of a solid body on a planet, it may be 
the upward jet and not the splash which is the domi- 
nant source of planetary material transferred into space 
|25|. Similarly, an oblique jet resulting from an oblique 
impact allows for an enhanced sidewards transport of ma- 
terial, as compared to the splash, (ii) The collapsing jet 
may contribute to the central peak often found in impact 
craters ji3:ll3- (iii) The downward jet will considerably 
change the layering of the sediments underneath a crater, 
as it provides a mechanism how surface material can be 
transported deep into the ground (see the supplementary 
material). In addition, a granular eruption will rearrange 
the sediment. Our suggested mechanism may shed new 
light on the sediment layering data found underneath the 
Chicxulub crater, which is a source of major controversy 

[illlilailp. 
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Similar peaks are observed in many craters of the terres- 
trial planets [lllTl 

The 3D simulations we did are too strongly affected by 
finite size effects. Nevertheless, also for these simulations 
a jet emerges. 

In our experiments we are in the region where the pres- 
sure increases linearly with depth |17| . 
In a geophysical context, exactly vertical impacts are of 
course very unlikely. 50% of all impacts on terrestrial 
planets occur with an angle between 30° and 60° |2^. 



